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Abstract 

We consider a 3-dimensional Dirac operator Ho with non-constant magnetic field of constant direction, 
perturbed by a sign-definite matrix-valued potential V decaying fast enough at infinity. Then we determine 
asymptotics, as the energy goes to +m and — m, of the spectral shift function for the pair (Ho, Ho + V). We 
obtain, as a by-product, a generalised version of Levinson's Theorem relating the eigenvalues asymptotics of 
Ho + V near +rn and — m to the scattering phase shift for the pair (Ho, Ho + V). 
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1 Introduction 



It is known [38] that the free Dirac Hamiltonian H m acting in the Hilbert space H := L 2 (]R 3 ; C 4 ) is unitarily 
equivalent to the operator h(P) ffi —h(P), where P := -iV and R 3 3 £ h-> := (£ 2 + m 2 ) 1/2 . For this 
reason, the set {±m} = /i[(V/i) _1 ({0})] of critical values of h plays an important role in spectral analysis 
and scattering theory for Dirac operators. For instance, one cannot prove at ±m the usual limiting absorption 
principle for operators H m + V, even with V a regular perturbation of H m , by using standard commutator 
methods. Both the statements and the proofs have to be modified (see e.g. [4, 19]). 

In this paper, we provide a new account on the spectral analysis of Dirac operators at the critical values by 
discussing the behaviour at ±m of the spectral shift function associated to sign-definite perturbations of Dirac 
operators with non-constant magnetic fields. Our work is closely related to [27] where G. D. Raikov treats a 
similar issue in the case of magnetic Pauli operators. It can also be considered as a complement of [33], where 
general properties of the spectrum of Dirac operators with variable magnetic fields of constant direction and 
matrix perturbations are determined. Other related results on the spectrum of 3-dimensional magnetic Dirac 
operators can be found in [2, 3, 5, 7, 11,13, 15, 16, 17, 20, 23, 34, 36, 37]. 

Let us describe the content of this paper. We consider a relativistic spin- 4 particle evolving in R 3 in pres- 
ence of a variable magnetic field of constant direction. By virtue of the Maxwell equations, we may assume 
with no loss of generality that the magnetic field has the form 

B(x\,x 2 ,x 3 ) = (0, 0, b(xi,X2))- 

The system is described in H by the Dirac operator 

H := aiUi + a 2 U 2 + a 3 P 3 + (3m, 

where (3 = ao, a±, a 2 , a 3 are the usual Dirac-Pauli matrices, m > is the mass of the particle and IX, := 
— idj — a,j are the generators of the magnetic translations with a vector potential 

a(xx, x 2 ,x 3 ) = (ai(x 1 ,x 2 ),a 2 (x 1 ,x 2 ),0) 

that satisfies B = dia 2 — d 2 ai. Since 03 = 0, we write P 3 = —id 3 instead of IT3. We assume that the function 
b : M 2 — > R is continuous (see Section 2 for details), so that Hq, defined on QJ°(IR 3 ; C 4 ), can be extended 
uniquely to a selfadjoint operator in H with domain T>(H^) . 

Then we consider a bounded positive multiplication operator V s C(R 3 ; ^h(C 4 )), where J?h(C 4 ) is the 
set of 4 x 4 hermitian matrices, and define the perturbed Hamiltonian H± := Hq ± V . Since V is bounded and 
symmetric, the operator H± is selfadjoint in TC and has domain T>{H) = V{Hn). We also assume that 
decays more rapidly than |.t|~ 3 as \x\ — > 00 and that 

(H± - z)- 3 - (H Q - z)- 3 G Sx(H) for each zel \ {<t(Hq) U <t(H±)}, (1.1) 

where S\(H) denotes the set of trace class operators in H. 

Under these assumptions, there exists a unique function £( ■ ; H±, Hq) € L 1 (R; (1 + | A|)~ 4 dA) such that 
the Lifshits-Krein trace formula 

Tr [f(H ± ) - f(H )] = [ d\ /'(A) £(A; H± , H ) (1.2) 

holds for each / G C^°(R) (see [39, Sec. 8.11]). The function £( • ; H ± ,H ) is called the spectral shift function 
for the pair (H±,Hq). It vanishes identically on R \ {o(Hq) U cr(H±)}, and can be related to the number of 
eigenvalues of H± in (—to, m) (see Remark 4.5). Morever, for almost every A G <r ac (Ho) the spectral shift 
function is related to the scattering matrix 5(A; H±, Hq) for the pair (H±,Hq) by the Birman-Krein formula 

detS(X;H ± ,H ) = c -^i(^H±Mo) 

After identification of £( • ; H±,Hq) with some representative of its equivalence class, our results are the 
following. In Proposition 4.4, we show that there exists a constant ( > defined in terms of b (cf. Proposition 
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2.1) such that £( • ; H± 7 Hq) is bounded on each compact subset of (—y/m 2 + £, y/ m 2 + £) \ {±to} and is 
continuous on ( — ^Jm 2 + ^/ m 2 + £) \ ({±m} U a p (iJ±)). In Theorem 6.5, we determine the asymptotic 
behaviour of £ ( A ; H± , Hq) as A — ► ±m, |A| < to, and in Theorem 6 . 1 4, we determine the asymptotic behaviour 
of £(A; H±, Ho) as A — > ±m, |A| > m. In both cases, one has ^(X;H±,Ho) — ► ±oo as A — ► =pm. The 
divergence of £(A; H±, Ho) near A = ±to scales as the number of eigenvalues near of certain Berezin- 
Toeplitz type operators. When V admits a power-like or exponential decay at infinity, or when it has a compact 
support, we give the first term of the asymptotic expansion of £(A, ; H±,Hq) near A = ±m (see Proposition 
6.10 and Corollary 6.17). In these cases, we show that the limits 

Z(m + e;H-,H ) , .. t(-m-e;H+,H ) 

lim — ) f and lim — ) f 

e\o £(to - e; H-,Ho) e\o £( - m + e; H+, H ) 

exist and are equal to positive constants depending on the decay rate of V at infinity (see Corollary 6.18 for a 
precise statement). This can be interpreted as a generalised version of Levinson's Theorem for the pair (H± , Ho) 
(see [21, 22] for usual versions of Levinson's Theorem for Dirac operators). The relation between the behaviour 
of the spectral shift function near A = +to and near A = —to is explained in Remark 6.15 by using the charge 
conjugation symmetry. 

These results are similar to the results of [27] (where Pauli operators with non-constant magnetic fields 
are considered) and [12] (where Schrodinger operators with constant magnetic field are considered). Part of the 
interest of this work relies on the fact that we were able to exhibit a non-trivial class of matrix potentials V 
satisfying (1.1) even though Ho is not a bounded perturbation of the free Dirac operator. We refer to Remark 
3.3 and Section 7 for a discussion of this issue. 

Let us fix the notations that are used in the paper. The norm and scalar product of 7i = L 2 (R 3 ; C 4 ) are 
denoted by j ■ || and ( • , ■ ). The symbol ® stands for the closed tensor product of Hilbert spaces and S P (H), 
p G [l,oo], denotes the p-th Schatten-von Neumann class of operators in H {Soo(Tt) is the set of compact 
operators in H). We denote by | • || p the corresponding operator norm. The variable x G R 3 is often written as 
X = (x±, X3), with x± G R 2 and x% G R. The symbol Qj, j = 1, 2, 3, denotes the multiplication operator by 
Xj in H, Q := (Qi, Q2, Q3), and Q± := (Qi, (^2)- Sometimes, when the context is unambiguous, we consider 
the operators Qj and Pj as operators in L 2 (R) instead of H without changing the notations. Given a selfadjoint 
operator A in a Hilbert space Q, the symbol E A ( • ) stands for the spectral measure of A. 

2 Unperturbed operator 

Throughout this paper we assume that the component b : M 2 — > M of the magnetic field B = (0, 0, b) belongs 
to the class of "admissible" magnetic fields defined in [27, Sec. 2.1]. Namely, we assume that b = bo + b, where 
60 > is a constant while the function b : R 2 — > R is such that the Poisson equation 

Aip = b 

admits a solution 92 : R 2 — > R, continuous and bounded together with its derivatives of order up to two. We 
also define ipo{x±) := i6o|xj_| 2 for each x± G R 2 and set tp := ipo + <P- Then we obtain a vector potential 
a = (ai,a 2 , 03) G C 1 (R 2 ; R 3 ) for the magnetic field B by putting 

a\ := d\(p, 0,2 '■= d2f and 0,3 := 0. 

(changing, if necessary, the gauge, we shall always assume that the vector potential a is of this form). We refer 
to [27] for further properties and examples of admissible magnetic fields. 

Since the vector potential a belongs to L ( g c (R 2 ; R 3 ), the magnetic Dirac operator 

H = ailli + a 2 n 2 + a 3 P 3 + [3m 

satisfies all the properties of [33, Sec. 2.1]. The operator Hq is essentially selfadjoint on C^°(R 3 ;C 4 ), with 
domain T>(Ho) C TL \j£ (R 3 ; C 4 ), the spectrum of Ho satisfies 

a(H ) = <7 a c(#o) = (-00, -m] U [to, 00), (2.1) 
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and we have the identity 

$l+l«i(P3 2 +m 2 ) \ 

n £7+ /cm i -\ tK^i E> 2 i ™ 2 n n 

(2.2) 



Hi 



( /T[«)l+l<»(P3 2 +m 2 ) 

ff+«)l + l«i(P|+m 2 ) 

H~ + (Ps+m 2 ) 

\ -frj«il+l®(Pf +rn 2 ) / 



with respect to the tensorial decomposition L 2 (]R 2 ) ® L 2 (]R) of L 2 (M 3 ). Here the operators H± are the compo- 
nents of the Pauli operator H± := © H± in L 2 (R 2 ; C 2 ) associated with the vector potential (eti, 02). 

We recall from [27, Sec. 2.2] that dimker(i/J) = 00, that dimker(iT^) = and that we have the 
following result. 

Proposition 2.1. Let b be an admissible magnetic field with bo > 0. Then = inf a(H±) is an isolated 
eigenvalue of infinite multiplicity. More precisely, we have 

dim ker ) = 00 and (0, ()cR \ v(Hj_), 

where 

C := 26 e- 2osc(t?) and osc(^) := sup £(xj.) - inf 

Finally, since (0, ()d \ o-(H ± ), we know from [33, Thm. 1.2.(d)] that the limits 

lim(Q 3 }-^ 2 (H -At ie^iQs)-^ 2 , v 3 > 1, (2.3) 

e\,0 



exist for each A G {—\Jm 2 + (, \Jm 2 + () \ {±to} (note that we use the usual notation ( • ) := + \ ■ | 2 ). 

3 Perturbed operator 

We consider now the perturbed operators H± = Ho ± V, where V = {Vjk} is the multiplication operator 
associated to the following matrix-valued function V. 

Assumption 3.1. The function V G C(R 3 ; ^(C 4 )) satisfies for each x = (x±,x 3 ) G M. 3 and each j,k G 
{!,... ,4} 

V{x) > and \ V jk (x)\ < Const. (xj_)^ l/± (x 3 )~ l/3 for some v± > 2 and v 3 > 1. (3.1) 
The potential V in Assumption 3.1 is short-range along x 3 . So we know from [33, Thm. 1.2] that 

(i) <7 e8S (i?±) = a css (H ) = (-00, -m] U [to, 00). 



(ii) The point spectrum of H± in ( — \J m 2 + (, \J to 2 + Q \ {±m} is composed of eigenvalues of finite 
multiplicity and with no accumulation point. 

(iii) H± has no singular continuous spectrum in ( — \J to 2 + £, •*/ to 2 + £) . In particular, and H± have a 
common spectral gap in (—to, to). 

Using the formula 

/>oo 

(A + A)- 7 = T( 7 )- 1 / dtp- 1 e-*t A+x > , A:V(A)^H, A>0, A, 7 >0, 
the diamagnetic inequality [1, Thm. 2.3], and the compactness criterion [9, Thm. 5.7.1], we find that 

l^l 1/2 (E £ <3 n l n ^ + TO2 ) _1/4 €5 00 [L 2 (]R 3 )]. 
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Since b is bounded this implies that 

\H,\-^V\H Q \-^ < |ffo|- 1/2 (E ilfc <4l^l)l^o|- 1/2 G SooCH). 

So |iJ r 1/2 V"|ff r 1/2 also belongs to Soo(W), since ^(W) is an hereditary C*-subalgebra of SSiJi) [24, 
Cor. 3.2.3]. One has in particular 

V^ 2 (\H Q \ + 1)- 1/2 e S^H). (3.2) 
The standard criterion [31, Thm. XI. 20] shows that 

|V^| 1/2 (- A + m 2 y G ^(R 3 )] if g G [2, oo) and iq > 3/2. 

This together with arguments as above implies that 

F 1/2 |# o r 7 G S q (H) if q > 2 is even and jq > 3. (3.3) 

So we have in particular that 

yi/2 E H G g 2 pq foj . any bounded borel set b CR. (3.4) 

In the sequel we shall need a more restrictive assumption on V. For this, we recall that there exists numbers 
\ {<j(Hq) U <t(H±)} since Hq and H± have a common spectral gap in (— m, m). We also set Rq(z) := 
(H - z)" 1 and R±{z) := (H± - for z G C \ tr(iJo) and z e C \ cr(H ± ), respectively. 

Assumption 3.2. The function V G C(]R 3 ; (C 4 )) satisfies for each x G ffi 3 anc/ eac/z j, fc G {1, . . . , 4} 

V(a;) > ant/ | Vjk(x)\ < Const. (x)~ v for some constant v > 3. (3.5) 

Furthermore, V is chosen such that 

R\{z) - Rl{z) G Si(H) foreachzeM\{a(Fo)U < T( J ff ± )}. (3.6) 

Note that (3.5) implies (3.1) if one takes v$ G (1, v — 2) and := v — v%. Note also that the choice of 
function A i— ► (A — z) -3 in the trace class condition (3.6) has been made for convenience. Many other choices 
would also guarantee the existence of the spectral shift function for the pair (H±, Hq) (see e.g. [39, Sec. 8.1 1]). 

Remark 3.3. Since the operator Hq is not a bounded perturbation of the free Dirac operator, we cannot apply 
the results of [40, Sec. 4] to prove the inclusion (3.6) under the condition (3.5). In general, one has to impose 
additional assumptions on V to get the result. For instance, ifV verifies (3.5), and 

(i) [V,a t ] = [V,a 2 ] = 0, 

(ii) for each i£l 3 and each j,k,£ G {1, . . . , 4}, one has \ (dg Vjk)(x)\ < Const. {x)~ <; for some <; > 3, 

(Hi) for each j, k, i G {1, . . . , 4}, one has (d e d 3 V jk ) G L°° (R 3 ), 

then (3.6) is satisfied. Furthermore, if V is scalar, then the same is true without assuming (Hi) (and (i) is 
trivially satisfied). The proof of these statements can be found in the appendix. Here, we only note that a matrix 
V G ^h(C ) satisfying (i) is necessarily of the form 

/ V! o V 3 _o_\ 

\/ — I v 2 V 3 I 

v I v 2 I ' 

\ V 3 V! / 

with Vx, V2 G R and V3 G C. 
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4 Spectral shift function 

In this section we recall some results due to A. Pushnitski on the representation of the spectral shift function for 
a pair of not semibounded self adjoint operators. 

Given a a Lebesgue measurable set B C R, we set p,{B) := ~ J B j^s, and note that p,(R) = 1. Further- 
more, if T = T* is a compact operator in a separable Hilbert space Q, we set 

n±(s;T) := rank^ ±T ((s, oo)) for s > 0. 

Then we have the following estimates. 

Lemma 4.1 (Lemma 2.1 of [26]). Let T x = T* G S 00(H) and T 2 = T 2 * G Si(H). Then one as for each 

S\,S2 > 

/ dn(t) n±(s 1 + s 2 ;T 1 +tT 2 )<n±(s 1 ;T 1 ) + — ||T a ||i. 

Jr 7rs 2 

For z£C \ (t(Hq), we define the usual weighted resolvent 

T(z) :=V l ' 2 (H Q -z)- l V 1 / 2 

and the corresponding real and imaginary parts 

A(z) := ReT(z) and := Im T(z). 

Then the next lemma is direct consequence of the inclusions (3.2)-(3.4) and [25, Prop. 4.4. (i)]. 

Lemma 4.2. Let V satisfy Assumption 3.1. Then, for almost every X G K, the limits A(X+iO) := lim e \^o A(X+ 
ie) and B(X + iO) := lim^o B(X + ie) > ex/sf ;n 84(H). 

Next theorem follows from the inclusions (3.2), (3.4), (3.6), from the equations (1.9), (8.1), (8.2) of [25], 
and from Theorem 8. 1 of [25]. 

Theorem 4.3. Let V satisfy Assumption 3.2. Then, for almost every X G ffi, £(A; H±, Hq) exists and is given by 
Z(X;H±,H Q ) = ± [ dfj,(t)n T (l;A(X + iO)+tB(X + iO)). (4.1) 



We know from (2.3) that A(X + iQ) and B(X + iO) exist in 8§(H) for each A G (-y/m 2 + (, y/m 2 + Q \ 
{±m}. In Propositions 5.2-5.3 and Corollary 5.5 below we show that in fact A(X+iO) G 84(H) and B(X+iO) G 
Si(H) for each A G (— \/m 2 + C, \Jm 2 + () \ {±m}. Hence, by Lemma 4.1, the r.h.s. of (4.1) will turn out 
to be well-defined for every A G ( — \/m 2 + £, y m 2 + Q \ {±to}. In the next proposition we state some 
regularity properties of the function 

(- Vm 2 + C, Vm 2 + \ {±m} bX^ £(A; fl±,.ffo) :=± / d M*) "t C 1 ? + i0 ) + tB ( X + M )) ■ 

The proof (which relies on Propositions 5.2-5.3, Lemma 5.4, Corollary 5.5 and the stability result [14, Thm. 3.12]) 
is similar to the one of [6, Sec. 4.2.1], 

Proposition 4.4. Let V satisfy Assumption 3.1. Then £( • ; H±, H ) is bounded on each compact subset of 
(— y/ m 2 + C, yj m 2 + C) \ {±m} and is continuous on (— \J m 2 + £, ^/ m 2 + £) \ ({±m} U c p (H±)) . 

In the sequel, we identify the functions £( • ; ifo) and £( • ; H±, Ho) since they are equal for almost 
every A G R due to Theorem 4.3 (see [35] for a study where the r.h.s. of (4.1) is directly treated as a definition 

af£(A;fl±,J7 )). 

Remark 4.5. In the interval (—to, to), i?o few «o spectrum and the spectrum of H± is purely discrete. Thus 
the spectral shift function £( • ; H±, Hq) can be related to the number of eigenvalues of H± as follows: for 
Ai , A2 G (—to, to) \ <r(H± ) with Ai < A2, we have ( see [25, Thm. 9.1 ]) 

Z(X 1 ;H±,H )-£(Xr,H±,H )=rar>kE H ±([X u X 2 )). 
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5 Decomposition of the weighted resolvent 

In this section we decompose the weighted resolvent 

T{z) = V 1 ' a {H -z)V 1 ' 2 , z€C\a(H ), 

into a sum T(z) = T d j v (z) + T^ oun ^{z), where T d i v (z) (respectively Tb oun d(z)) corresponds to the diverging 
(respectively non-diverging) part of T(z) as z — > ±m. Then we estimate the behaviour, in suitable Schatten 
norms, of each term as z — > ±m. We refer to [12, Sec. 4] and [27, Sec. 4.2] for similar approaches in the case 
of the Schrodinger and Pauli operators. 

Let a and a* be the closures in L 2 (M 2 ) of the operators given by 

a<p := (III — «n 2 )( ( 9 and a*<p := (Hi + illa)^, 

for ip e C§°(K 2 ). Then one has (see [38, Sec. 5.5.2] and [28, Sec. 5]) 

/ m 1®P 3 a®l \ 
tt I m a*®l -1®P 3 I /c n 

■"0 IfflP, a (9.1 _m nil \ J -LJ 



with 



Now, let 



l(g>P 3 a®l —to 
\ a*®l -1®P 3 -m / 

ker(a*) = ker(aa*) = ker(ifj) C L 2 (M 2 ). (5.2) 

P 

P 


be the orthogonal projection onto the union of the eigenspaces of Ho corresponding to the values A = ±m. 
Since P = p (g> 1 is the orthogonal projection onto ker(ifj) <8> L 2 (R), the equations (5.1) and (5.2) imply that 
Hq and P commute: 

H^P = PH^. (5.3) 
In fact, by using (2.2) and (5.1), one gets for each z E C \ <r(i?o) the equalities 

(H^-z)- 1 ? 

= (H a + z)(H 2 -z 2 y 1 P 

( (z+to) \ / 1 \ 

= [p®R(z 2 - m 2 )}( § o o o +[p®P 3 ^ 2 -m 2 )] ? 8 8 8 . 

V o o o 0/ Voooo/ 

where R{z) := (P$ - z)~ ,z e C\[0, oo), is the resolvent of P| in L 2 (M). This allows us to decompose T(z) 
as T(z) = T div (z) + T bound (z), with 



/ (z+m) 0\ 

T d Uz):=V^[p®R(z 2 -m*)}( g go g ^, 



0. 



T bound (z) := F 1 / 2 [p (X) P^R(z 2 - m 2 )] f ° g ° V 1 ' 2 + V 1 ^ 2 (H - z^P^V 1 / 2 (P x := 1 - P). 

Voooo/ 

One may note that this decomposition of T(z) differs slightly from the simpler decomposition 

T(z) = V^iHo - z)PV 1 ' 2 + V l ' 2 (H - z)P^V 1 ' 2 , 

since the first term in Tb OU nd(z) is associated to the projection P and not the projection P^-. This choice is 
motivated by the will of distinguishing clearly the contribution T d j v (z), that diverge as z — ► ±m, from the 
contribution Tt >oun d(z), that stays bounded asz^ ±m. 
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For A £ R \ {0}, we can define the boundary value i?(A) of the resolvent R(z) as the operator with 
convolution kernel r\( ■ ), where 



r\(x 3 ) := 



i7=sr ifA<0 > 



ifA>0, 



for each £3 £ R. So, we can extend the definition of Tdi v ( ■ ) to the values A G 1\ {±to}: 

((A+m) \ 
0000/ 

In the following proposition, we show that the trace norm of Tdi v (z) is continuous in C+ := {z e C 
\m(z) > 0} outside the points z = ±m, where it may diverge as \z =F m| -1 / 2 . The proof of the proposition 
relies on a technical result that we now recall. 

Lemma 5.1 (Lemma 2.4 of [27]). Let U £ L 9 (R 2 ), q £ [1, 00), and assume that b is an admissible magnetic 
field. Then pUp £ S q [L 2 (R 2 )}, and 

\\r>TJr>\\ q < — p 2osc (^) Hr/11 9 

The symbol y + denotes the postive part of y £ R. 
Proposition 5.2. Let V satisfy Assumption 3.1. Then the operator-valued function 

C+ \ {±m} 3 z ^ T div (z) £ S^H) 
is well-defined and continuous. Moreover, we have for each A £ R \ {±m} the bound 

ll^iv(A)|U < Const. (||±m| 1/2 + l^ 2 ) (l + (A 2 - m»#*). 

Proof. We have for each z £ C \ ct(Hq) the identity 

((z+m) 0\ 
0(A) 8 )M, 
0000/ 

where 

M := F 1/2 (Q±)^ /2 (Q 3 P /2 , (5.4) 
G := (Q±)- /J - /a p{Q±)- Vj - /2 , (5.5) 
J, := (Q3K 3/2 i?M(Q3)- l ' 3/2 . 
The operator M is bounded due to Assumption 3.1. So 

||Tdiv(z)||i < Const. (\z + m\ + \z-m\) \\G\\i\\ J z 2- m * ||i. 

But we know from Lemma 5.1 that ||G||i < Const., and from [6, Sec. 4.1] that the operator-valued function 
^+ \ {0} 3 z 1-^ J z is continuous in the trace norm and admits the bound 

||J A ||i < Const. (1 + A+ /4 )|Ar 1/2 , A e R\{0}. 

It follows that 

\\T div (z)h < Const. (|A±i|| 1/2 + \m 1/2 ) (1 + ~ - 2 )f) 
for each A G R\ {±m}. □ 
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In the following proposition, we show that the function z i— > T^ oun ^(z) £ Si(7i) is continuous in C \ 
{ (-00, -\Jm 2 + C] U [^m 2 + C, 00)} . The symbols ff* stand for the operators if* := iff (g> 1 + 1 ® P 3 2 
acting in L 2 (R 3 ). 

Proposition 5.3. Let V satisfy Assumption 3.1. Then the operator-valued function 



\ {(-oc,-Vm 2 +C] U [v/m 2 + C,(^)} 9 z ^ T bound (z) g fi*(W) 



is well-defined and continuous. Moreover, we have for each A G (— yjrn?~- f~£, ^/ m 2 + £) f/ze bound 

||T bound (A)|| 4 < Const. (|A| + A 2 ) (l + { *~™^$ + ) + Const. (5.6) 
Proof. One has the identity 

(Ho - z) -1 = #0"' + z(l + zHv 1 ) {H 2 - z 2 )- 1 
for each zeC \ a(Ho). Thus the operator T boU nd(2) can be written as 

T bound (z) = M(G ®S Z ) (VuI)m + V^H^V 1 ' 2 + zV^ 2 (l + zH^)(H 2 * 2 )" W 1 / 2 



(5.7) 

^T^+Tz+niz), 
with M and G given by (5.4)-(5.5), and 

S z := (Q 3 )" l ' 3 / 2 P 3 i?(z 2 -m 2 )(Q 3 )-^/ 2 . 

The integral kernel of S 2 is 

|(x3)-^/ 2 jg^?fe i ^ 5 =^l-3-4l(4)-3/2 ) (58) 



with the branch of \J z 2 — m 2 chosen so that Im \/ z 2 — m 2 > 0. So S z extends to an element of S^L^R)] for 
each z E C, with ||5 Z ||2 < Const. Since M is bounded and ||G||i < Const., this implies that 

||Ti(z)|| 2 < Const. ||M|| 2 ||G||i||5,|| 2 < Const. (5.9) 

for each z € C. One also has 

||r 2 || 4 < Const. (5.10) 

due to (3.3). So, it only remains to bound the term T^{z). 

Letz e C\{(-oo,- v /?ti 2 + C]U [yjm 2 + C,oo)} and P^ := 1- P. Then (H~ + m 2 - z 2 )^ 1 P 1 - and 
(H+ + m 2 - z 2 ) -1 belong to J?[L 2 (1R 3 )], and we have 



[H~ +m 2 - z 2 ) 1 P 1 - = P- L (H~ +m 2 - z 2 ) \ 



Thus 



(H 2 z 2 )~ W 1 / 2 = (H 2 z 2 )- 1 ("HID V" 2 



0001 

' P^(H~ +m A -z 2 )~ 1 

(H + +m 2 -z 2 )- 1 \ v l/2 

P ± (H-+m 2 -z 2 )- 1 
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and 

\\{H 2 - 2 2 )- W^HJ < 2\\M\\ 2 {\\P^(H- + m 2 - z 2 )~ lM 4l + \\{ H+ + ™ 2 - z 2 ) -1 M 2 \\\) , 

where M 2 := {Qi_)~ Vx/2 (Q 3 )~ V3/2 . But, we know from the proof of [27, Prop. 4.4] that 

H-P-^//"" +m 2 - z 2 )- 1 M 2 \\ 2 < Const. C(z) and ||(iJ+ + m 2 - z 2 )- 1 M 2 \\ 2 < Const. C(z), 
where 



y + 1 

iup 

It follows that 



CO) := sup | — ■ ^ 2j 

[Coo) |y + m - ^ I 



||r 3 («)|| a < Const.||zF 1 / 2 (l + 2 J ff - 1 )|| p/||C*(z) < Const. (|z| + |z| 2 ) C(z). (5.11) 

The claim follows then by putting together (5.9), (5.10), and (5.11). □ 

In the next lemma we give some results on the imaginary part of the operator S z in L 2 (R) appearing in the 
proof of Proposition 5.3 

S z = (Q 3 )-" 3 / 2 P 3 R(z 2 - m 2 )(Q 3 )-»^ 2 , zeC\ a(H ), u 3 > 1. 

Lemma 5.4. (a) One has \mS\ = Qfor each X £ (— m, m). 

(b) Let p > 1 be an integer. Then one has for each A £ R with |A| > m 

II lm5 A || p < C p , 
where C p is a constant independent of X. Furthermore 

lim ||lmS A ||p = 0. 

A — >±m, | A| >m 

Proof, (a) This is a direct consequence of the spectral theorem. 

(b) Let A £ R, |A| > m. Then one shows by using (5.8) that Im S\ is equal to the rank two operator 

lm 5 A = (v x , ■)u\+ (ma, ■ ) v\, 

with 

u\(x 3 ) := (xsy^sm (x 3 VX 2 -m 2 ) and v x (x 3 ) := -%(x 3 )-*' 2 cos (x 3 VX 2 -m 2 ). 
Since (v\, u\) = 0, this implies that 

| \ m S x \ p = \\u x \\ p (v X: -)vx + \\vx\\ p (ux, -)u x . 

Thus 

|| \mS x \\* = Tr (| Im S x \ p ) = KH P \\v x \\ 2 + |MI P Kl| 2 - 
This, together with the equality 

lim ||ua|| = 0, 

A — >±m, | A|>m 

implies the claim. □ 
In the next corollary we combine some of the results of Propositions 5.2, 5.3 and Lemma 5.4. 
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Corollary 5.5. Let V satisfy Assumption 3.1. Then the identity 

^(A + ^o) = ^ div (A) + ^ bound (A) (5.12) 

holds for each A £ (— \Jm 2 + £, \J m 2 + \ {±m}, anc/ the estimate 

|| lmT bound (A)|| p < Const. || lmS A j| p (5.13) 



holds for each integer p > 1 an each A £ (— \/m 2 + C, yV+C). /« particular, we have 



lim lmT bound (A) = 0, (5.14) 



due to Lemma 5.4. 



Proof. The first identity follows from Propositions 5.2 and 5.3. Let A £ ( — y^m 2 + £, m 2 + £). Using (5.7) 
and the commutation rule (5.3) one obtains that 

\mT bound (\) = M(G®\mS x ) (Vul) M, 

\0 0/ 

withM and G defined by (5.4)-(5.5). Since M is bounded and ||G||i < Const., this implies (5.13). □ 

6 Proof of the main results 

We begin this section by showing that the value of £(A; H, H± ) as A — > ±m is bounded from below and from 
above by expressions involving only the term T d i v (A) of the decomposition T(A + iO) = T d i v (A) + Tbound(A). 
Then we consider separately the limits A — > ±m with | A| < m and the limits A — > ±m with |A| > m. 

We start by recalling two standard properties of the counting functions n± . Given two compact operators 
T\ = T* and T2 = T 2 * in a separable Hilbert space Q, we have the Weyl inequalities 

n±(si + s a ;Ti+T a )<n±(*i;Ti)+n±(s a ;T 2 ) for each Sl , s 2 > 0. (6.1) 

Moreover, if T = T* belongs to S P (Q) for some p £ [1, 00), then 

n ± {s;T) < s~ p \\T\\ p p foreachs>0. (6.2) 

Proposition 6.1. Let V satisfy Assumption 3.2. Then the estimates 

dfi(t) n± (1 + e; ReT div (A) + t Im T div (A)) + 0(1) 
<T£(A;ff T ,# ) 

< / d/i(t)n ± (l-£;ReT div (A)+tlmT div (A)) + 0(1) 
hold as A — > ±m for each e £ (0, 1). 

Proof. Using (5.12), the Weyl inequalities (6.1), and Lemma 4. 1 we get 

/ d^(t)n ± (l + e;ReT div (A)+tlmT div (A)) - n T (e/2; ReT bound (A)) - — 1| Im T bound (A) || 

JR ^ 

dfi(t) n± (1; A(A + iO) + £B(A + iO)) 
< / dM(t)n ± (l-e;ReT div (A)+tlmT div (A))+n ± (e/2;Rer bound (A))+— ||lmT bound (A)|| (6.3) 



J 
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Due to (6.2), we have 



7i ± (e/2;ReT bound (A)) < 16£- 4 ||T bound (A)|| 4 , 



which combined with (5.6) gives 

n ± (e/2;ReT bound (A)) = 0(1) as A — > ±m. 
Moreover, we know from (5.14) that 

lim II lmT b0 und(A)|| = 0. 

So the claim follows from the estimates (6.3) and Formula (4.1) □ 

6.1 The case |A| < m 

In this section we prove asymptotic estimates for £(A; H, H±) as A — * ±m with |A| < to. We start with a 
corollary of Proposition 6.1, which follows from the fact that lmTdi v (A) = and ReTdi v (A) = Xdiv(A) for 
A G (—m, to). 

Corollary 6.2. Let V satisfy Assumption 3.2. Then the estimates 

n± (1 + s; Tdiv(A)) + 0(1) < t£(A; H t ,H ) < n± (l - e; T div (A)) + 0(1) 
/10/c/ ai A — > ±m, | A| < to, /or eac/i e G (0, 1). 

Define the bounded operators K± : H — > L 2 (M 2 ; C 4 ) by 

(i^)(*±) := / dx' ± dx' 3 p(x±,x'±.) ( 8 8 11) V 1 '*^, 4M*x, 4), 
Jm 3 Voooo/ 

(i^)(s±):= / dx^K^^fHHV^^'^M^^s), 

where p( • , • ) is the integral kernel of the projection p. One shows easily that K± : L 2 (R 2 ; C 4 ) — > 7i are given 
by 

(X; V) (XX , *3 ) = v l ' 2 (X X , *a) f 8 8 8 8 ) (W>) (*X ) , 

Voooo/ 

(A-^)(xx,x 3 ) = y 1 / 2 ( a ;x,x 3 )f88?n (WOK), 

Voooo/ 

and that 

O+W := i(^l) 1/2 ^;^ + ^nd 0-(A) := -§(^) V2 ^- 
belong to S^W) for each A G (— to, to). 

In the next proposition we show that the functions n± (■ ; T d j v (A)) as A — > ±m, |A| < 7n, can be bounded, 
up to 0(1) terms, from below and from above by expressions involving 0± (A). 

Proposition 6.3. Let V satisfy Assumption 3.2. Then the estimates 

n+((l + e)«;0+(A)) +0(1) < n+(s;T div (A)) < n+((l - e)s; 0+(A)) + 0(1), (6.4) 
0(1) <n_(s;T div (A)) < 0(1), (6.5) 

hold as \ Z' m, for each e G (0,1) a«t/ s > 0, anc/ the estimates 

0(1) <n+(s;T div (A)) < 0(1), (6.6) 
n_ ((1 + e)s; 0_ (A)) + 0(1) < n_ (s; T div (A)) < n_ ((1 - e)«; 0_ (A)) +0(1), (6.7) 

hold as A \ —m, for each e G (0, 1) anrf s > 0. 
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Proof. We only give the proof of (6.4)-(6.5), since the proof of (6.6)-(6.7) is similar. In point (i) below we 
show that the difference T d iv(A) — 0+(A) can be approximated in norm, as A / m, by a compact operator 
independent of A. Then we prove (6.4)-(6.5) in point (ii) by using this result, 
(i) Let A G ( — T7i, m) and take v' G (3, v). A direct calculation shows that 

~ / (A+m) °\ ~ 

T div (A) - 0+(A) = M{G V . V , ® g g o g M + 0_(A), (6.8) 



(A-m) 
0, 



where : L 2 (M) -> L 2 (R) is given by 



K> Vjfe) := -(a*) 7 / dr 3 smh( Ux 3 ) ' 

Jm. \m — A- v * 7 



v<(4) 



3^ 

and 

M := ^ 1 / 2 (g ± )^- !/ ')/ 2 <Q 3 )^'/ 2 , (6.9) 
:= (Q±)- {v -' / ' )/2 p(Q±)- i ^' )/2 . (6.10) 

The operator M is bounded due to Assumption 3.2, G„_„< is compact in L 2 (R 2 ; C 4 ) due to Lemma 5.1, and 
0_ (A) satisfies 

lim ||O_(A)|| a = 0. (6.11) 

A — >m, I A| <m 



Define 



/ (m±m) \ 

r ± :=M(G„_^jH) § g.^g M , (6.12) 



0, 



with : L 2 (R) -> L 2 (R) given by 

(4 m V)(^3) :=-i(^ 3 )- v ' /2 / dx' 3 \x 3 ~x' 3 \(x' 3 )^'^(x' 3 ). 



Since // > 3, jjy belongs to 52[L 2 (R)], and T± is compact in Ji. Moreover, by using Lebesgue's dominated 
convergence theorem, one shows that 

lim ||4 m) - 4^11" = 0. 
A— >±m, |A|<m 2 

This, together with (6.8), (6.1 1) and (6.12), implies that 

lim ||T div (A) - 0+(A) -2+11= 0. (6.13) 

A /^m 

(ii) Take A G (—771, m), £ G (0, 1), and s > 0. Using the Weyl inequalities (6.1) we get 

n± ((1 + e)s; 0+(A)) - n T (es; T div (A) - 0+(A)) < n± (s; T div (A)) 

< n± ((1 - e)s; 0+(A)) + n± (es; T div (A) - 0+(A)) . 

Now we have n_ (i; + (A)) = for each t > and A G (—m, to), since 0+(A) is a positive operator. So, to 
prove (6.4)-(6.5), it is sufficient to show that n± (es; T div (A) - 0+(A)) = 0(1) as A / m, for each e G (0, 1) 
and s > 0. Let t > be fixed. Then we know from (6.13) that we can chose A + G (—971, to), close enough to 
m, so that ||T<)j V (A+) — 0+(A + ) — T+|| < t/2. Thus, using again the Weyl inequalities, we get 

n±(t;T div (A + ) - 0+(A+)) < n± (i/2; T div (A+) - 0+(A+) - T+) + ri± (t/2; T+) = n± (t/2; T+) . 

Since the r.h.s. is independent of A + we have shown that n±(i;T d i v (A) — 0+(A)) = 0(1) as A / m. This 
concludes the proof of (6.4)-(6.5). □ 
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We show now that the counting functions n± ( ■ ; 0± (A)) in Proposition 6.3 can be rewritten in terms of 
Berezin-Toeplitz type operators. Define for each A G (—m,m) 

:= h{£3t) 1/a pW + P and W _(A) := -l(^) 1/2 p W.p, 
where the functions W± : M 2 — > R are given by 

W+(o:j.) := / dx 3 Vii(x_L,a;3) and W_(:zi_l) := / dx 3 V 33 (xj_, x 3 ). (6.14) 

Under the condition (3.5) one has 

< W±{x 1 _) < Const. (x±_)- u+1 for all x± G K 2 , 
andw±(A) G Si[L 2 (]R 2 )] if 1/ satisfies Assumption 3.1 (see Lemma 5.1). Moreover, one has the following. 
Proposition 6.4. Lef V satisfy Assumption 3.1. Then we have for each A G (—to, to) and s > 

n±(s;0±(A)) =n±(«;w±(A)). (6.15) 
Proof. Given s > and two separable Hilbert spaces 7^1,7^2, one has 

n± (s; = n± (s; SB*) (6.16) 

for any G 9B^H\ , 7i 2 ) such that B* B G Soo (?i 1 ) . Moreover, one can easily check that 

/ 1 o\ /o o\ 

A+A + = (8880 P w +P and K - K - =0010 PW-P. 
Voooo/ Voooo/ 

Thus 

n + (s;0 + (X)) = n + (.s; i(^±A) 1/2 (j jj jj jj) P^+p) = n + («; ± (S^/'pW+p) = n + (s;o, + (A)). 

The proof of the second equality in (6.15) is similar. □ 

The next theorem is direct consequence of Corollary 6.2 and Propositions 6.3-6.4. 

Theorem 6.5. Let V satisfy Assumption 3.2. Then one has for each e G (0, 1) 

0(1) <£(\\H + ,H ) <0(1) (6.17) 
-n+(l-e;w+(A)) +0(1) <£(\;H-,H Q ) < -n+(l + e;w+(A)) + 0(1) (6.18) 

as A /* to, ana" 

n_(l+e;w_(A))+O(l)<f(A;fr+,F )<n_(l-e;w_(A))+O(l) (6.19) 
0(1) < £(\;H-,H ) < 0(1) (6.20) 

as A \ —to. 

Remark 6.6. 77ze inequalities (6.17) together with Remark (4.5) imply that the eigenvalues of Hq + 1/ in 
(—to, to) near +to (if any) do not accumulate at +m. (9« fne of/ier /lana" the inequalities (6.18) tell us that the 
number of eigenvalues of Hq — V in (—to, to) near A = +m scales, up to 0(1) terms, as 



i+(s;uj + (X)) =r&nk Ep w +p(( 



*<^) 1/2 ,- 



vw'fn s ~ 2. Accordingly, the problem of counting the number of eigenvalues of Hq — V in (— to, to) near 
+to reduces to the problem of counting the number of eigenvalues of the positive Berezin-Toeplitz type operator 
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pW+p nearO. The inequalities (6.19)-(6.20) lead to similar conclusions on the number of eigenvalues of Hq±V 
in (— m, m) near — m. 

One can compare these results with the results of [8] and [19] on the finiteness in (— m, m) of the discrete 
spectrum of the Dirac operator perturbed by a matrix potential Q = {Qjk(x)}^ ^ =1 . In Corollary 2.2 of [8], the 
author shows that the spectrum in (—to, to) of the Dirac operator perturbed by Q is finite if the 2x2 diagonal 
blocks of Q are of order 0(\ i x\~ 2 ~ & } and the anti-diagonal blocks are of order for some 5 > as 

\x\ — > oo. In Corollary 2.1 of [19], the authors show that the Dirac operator perturbed by jQ, with \~/\ small 
enough and 

\Qjk{x)\ < (x)- 2 , j,k€ {1,2,3,4}, 

does not have any point spectrum. Therefore, in our case where Q = — a\a\ — ct2d2 + V, we would not have 
had any accumulation of eigenvalues in (— m,m) if we would have imposed such decay assumptions on the 
magnetic part —a\a\ — 0120-1 of the perturbation. 

As seen in Theorem 6.5 the behaviour of the function £( • ; H± , Hq) in (— m, m) depends on the distribution 
of eigenvalues of the trace class operator pW T p. In our next proposition we shall exhibit different types of 
behaviours depending on the choice of the functions Vu and V33 appearing in W±. For that purpose, we first 
have to recall some technical results taken from [27], [29] and [30]. 

In the first lemma, an integrated density of states (IDS) for the operator iJj in L 2 (R 2 ) is defined as follows 

(see e.g. [10, 18]): Let xt,x ± be the characteristic function of the square x± + ( — -j, -j) 2 , with x± G M 2 and 
T > 0. Then a non-increasing function g : [0, 00) — > R is called IDS for the operator i?J if for each x± <G R 2 
it satisfies 

£?(A) = lim T- 2 Tr[ X T, x AQ±)E H ^{(-™,\))xT, x AQ±-)\ 

1 — >oo 

for each point A £ R of continuity of g. 

Lemma 6.7 (Lemma 3.3 of [27]). Let U £ C 1 (R 2 ) satisfy 

< U(xj_) < Const. (x±}- a and \(VU)(xj_)\ < Const. (x 1 _)- a ~ 1 
for all x £ R 2 and some a > 0. Assume moreover that 

• U(x±) = + o(l)) as \x±\ — > 00, where u is a continuous function on S 1 which does not 
vanish identically, 

• b is an admissible magnetic field, 

• there exists an IDS gbfor the operator H^. 
Then we have 

n+(s; P Up) =^|{nel 2 | U{x x ) > fl}|(l + o(l)) = * Q (s; u, b Q ) (l + o(l)) as s \ 0, 
where \ ■ | denotes the Lebesgue measure, and 

^ a (s;u,b ) ■= S T" 60 / di?w(i?) 2/Q , s>0. (6.21) 
4tt J s i 

Lemma 6.8 (Lemma 3.4 of [27]). Let < U £ L°°(R 2 ). AiiMme that 

In (l/(ix)) = -r/|x_L| 2/3 (l + o(l)) as 00, 
/or some r], (3 > 0. Lef 6 oe an admissible magnetic field. Then we have 

n + (s\pUp) =$ / j(s,?7,&o) (l + o(l)) as s \ 0, 
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where 

r^viin^ 1/0 6(0,1), 

*f,{8,r,M) ■= \ m(x+2 V / bo ) \Hs)\ '7 = 1, sefO.c- 1 ). (6.22) 

[^(lnlln^D-'lMs)! ifP>l, 

Lemma 6.9 (Lemma 3.5 of [27]). Let < U € L°°(]R 2 ). Assume that the support ofU is compact, and that 
there exists a constant C > such that U > C on an open non-empty subset of M 2 . Let b be an admissible 
magnetic field. Then we have 

n+(s;pUp) = $ DO (s)(l + o(l)) as s \ 0, 

vv/tere 

*„,(«) :=(]n|M«)l) _1 |M«)l. sefO.c- 1 ). (6.23) 

Combining Theorem 6.5 with Lemmas 6.7-6.9 we obtain the behaviour of ^(X;H±,Ho) as |A| — » to, 
|A| < m, when the functions W± admit a power-like or exponential decay at infinity, or when they have a 
compact support. 

Proposition 6.10. Let V satisfy Assumption 3.2. 

(a) Assume that the hypotheses of Lemma 6.7 hold with U± = W± and a = V — 1. Then we have 

a\;H-,H Q ) = -*,- 1 (2(2i T A) 1/2 ; U+ ,&o)(l + (l)) as X / to, 

and 

C(A ;j ff + ,i/ ) = ^-i(2(^) 1/2 ; U -,6 )(l + o (l)) as A \ — m, 
vw'f/i ^^-l g;Ven by Equation (6.21). 

(b) Assume that the hypotheses of Lemma 6.8 hold with U± = W±. Then we have 

e(A;i/_,ifo) = -$ 3+ (2(^) 1/2 ;? 7+,&o)(l + o(l)) as X / m, 

and 

C(A; J ff+,i7 ) = $ /3 _(2(2i±A) 1/2 ;r ? _,&o) (l + o(l)) as X \ —to, 
with j3± € (0, oo) and $>/3 ± given by Equation (6.22). 

(c) Assume that the hypotheses of Lemma 6.9 hold with U± = W±. Then we have 

e(A;i/_,H ) = -$oo(2(^) 1/2 )(l + O (l)) as X / to, 

and 

C(A;i7 + ,i/ ) = $ co (2(f±A) 1/2 )(l + O (l)) as X \ — m, 

with $oo given by Equation (6.23). 

The estimates of Proposition 6.10 are similar to the ones of [27, Cor. 3.6], where the corresponding situation 
for magnetic Pauli operators is considered. 
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6.2 The case |A| > m 

In this section we prove asymptotic estimates for £(A; H, H±) as A — > ±m, when |A| > m. We start by showing 
an estimate for n± (s; ReT d j v (A)). 

Proposition 6.11. Let V satisfy Assumption 3.2. Then the estimates 

n±(s;ReT div (A)) = 0(1) as A — > ±m, |A| > m, 

hold for each s > 0. 

Proof. Take A € K with |A| > m, and let i^' e (3, f ). Then we have 

ReT div (A) = M(G V - V i ® i#>) § g (A _° m) o M, 

\ 0/ 

with M and as in (6.9)-(6.10), and 

B$ := (Q 3 )- 1 ' 72 Rei?(A 2 - m 2 )^}"" 72 . 
By using Lebesgue's dominated convergence theorem, one shows that 

Urn II ReT div (A) - T ± \\ = 0, 

A — ► im, | A| >m 

with T± as in (6.12). So the claim can be proved as in point (ii) of the proof of Proposition 6.3. □ 

The next result follows from applying Propositions 6. 1 and 6. 1 1, the Weyl inequalities (6. 1) and the iden- 
tities [12, Sec. 5.4] 

/ d^(t)n ± (s;tT) = 7 r- 1 Trarctan(s- 1 T), s > 0, (6.24) 

JR 

where T £ Si(H), T = T* > 0. We also use the fact that sgn(A) Im T d iv(A) is apositive operator if |A| > m. 
Corollary 6.12. Let V satisfy Assumption 3.2. Then the estimates 

ir- 1 Trarctan [(1 + e)" 1 sgn(A) Im T div (A)] + 0(1) 

< T£(A;ff T ,ffo) 

< ir- 1 Trarctan [(1 - e)" 1 sgn(A) Im T div (A)] + 0(1) 

hold as A — > ±m, |A| > to, for each e £ (0, 1). 

As in the case | A | < to, we introduce auxiliary operators in order to express the lower and upper bounds for 
=F£(A; H t , H ) in terms of Berezin-Toeplitz type operators. For A G M with | A| > m, we define the operators 

K lt x,K a>x : L 2 (M 2 ;C 4 )by 



y/\X+m\ 



(K ltX <p){x ± ) := j dx' ± dx' s p(x ± ,x' ± ) C o S (x' 3 V>^?) [ ° o ^_ o j ^(^,4)^,4), 



0, 



y/\X+m\ 



(K\ W )(x ± ) := / d^d4K^,^)sin(4v / A^^) ( ° ° o ] V 1 / 2 ^, x' 3 Mx' ± , x' 3 ). 



0, 
2fn2.iri\ 



Direct calculations show that the adjoint operators ^, K% \ : L (R ; C ) — > W are given by 



^/ 1 A+m | 



(KlMx ± ,x 3 ) = cos(x 3 \/A 2 -m*)vV 2 {x±,x 3 ) \ ° ° ° ) (^)(z±), 



, 



yfA+mj 



(KlMx ± ,x 3 )= S m{x 3 y/>?-Tn 2 )V 1 / 2 (x ± ,x 3 ) ° ° °- ° (pi&Xsj.), 

y|A— m| 1 



0, 
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and that 

sgn(A) lmT div (A) = — — L= (A* A A liA + K* X K 2<X ). 
2v \ l — m z 

This last equation can be written more compactly as 



sgn(A) Im T div (A) = A* A A (6.25) 

2\J \ z — m z 



if we use the operator 



K X :H^L 2 (R 2 ;C 8 ), K x ip :-- 

with adjoint 



K* x : L 2 (M 2 ; C 8 ) -> W, A'J (Jj - A^i + # 2 * A ^2- 

For the next proposition we also need to introduce for each A £ R with |A| > m the positive operator 
0(A) : L 2 (R 2 ;C 8 ) -> L 2 (R 2 ; C 8 ) defined by 

A direct calculation shows that 



#i,A<A 



A A A* = p 
where 



M 2 , A M 3 , A 



Mi jA (a;j_) := / da; 3 cos 2 (x 3 \/A 2 - m 2 ) 



|A+m|V r i 1 (i J .,x 3 ) VA^m^VisCxx.is) O" 



VA^-m^Vbi^j.,^) \X-m\V 33 (x ± ,x 3 ) 

0/ 



/ |A+m|Vu(xj_,a;3) y/X*-m- i V 13 (x ± ,x 3 ) 

M 2)A (a;_ L ) = / dxa sin (x 3 \/A 2 - to 2 ) cos (x- 3 \/A 2 - m 2 ) ' 



VA^m^siln^s) |A-m|y 33 (2;x,X3) I ' 
0/ 



M: 



|A+m|Vu(a:x,a!3) V A^ -m J V 13 (x ± ,x 3 ) 

m It/ C \ n 11 n/° f \ n 

Its \ VX 2 -m 2 V 31 (x±,x 3 ) | A— m| V33 (li ,a; 3 ) 

0/ 



iejJ = / dx 3 sin 2 (x 3 \/ A 2 — m 2 ) 
This implies that 



||0(A)|| a < (j^^WpW+pl + ($3%) 1/2 \\pW-p\\ v 

and thus 0(A) G Si [L 2 (R 2 ; C 8 )] if V satisfies Assumption 3.1. 

Next Proposition is a direct consequence of Equations (6.16) and (6.25). 

Proposition 6.13. Let V satisfy Assumption 3.1. Then we have for each A G R with |A| > to and each s > 

n ± (s; sgn(A) Im T div (A)) = n± (s; 0(A)) . 

/« particular, it follows by Equation (6.24) that 

Trarctan (s _1 sgn(A) Im Tdiv(A)) = Trarctan (s _1 0(A)). (6.26) 

The combination of Corollary 6.12 and Equation (6.26) gives the following. 

Theorem 6.14. Let V satisfy Assumption 3.2. Then one has for each e € (0, 1) 

iTT" 1 Trarctan [(l±e)- 1 0(A)] +0(1) < ^(X;H ± ,H ) < +7T" 1 Trarctan [(1 =p e)~ 1 0(A)] +<D(T) 

as A — > ±m, | A| > m. 
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Remark 6.15. The fact that the operators u± (A) and f2(A) in Theorems 6.5 and 6.14 depend in a distinguished 
way on the components V\\ and V33 ofV is due to our initial assumption 60 > 0. Indeed, this choice implies 
that ker(TJj) is non trivial, whereas ker(TJ^) = {0}. This lead us to introduce in Section 5 the projection 
P = diag(P, 0, P, 0), which put into light the priviledged role of the components V\\ and V33 ofV. 

The variation of £(A; H± , Hq) under the change A 1— ► — A can be explained using the antinunitary trans- 
formation of charge conjugation [38, Sec. 1.4.6] 

where Uc '■= ifict^. Indeed, if we write H(a, ±V) and H q{o) for H± and Hq, then a direct calculation using 
the Lifshits-Krein trace formula (1.2) shows that 

Cfl"(o,±V)C _1 = -H(-a,TU c VUc), 

which entails 

f (A; H(S, ±V),H (Zj) = -£( - A; H(-S, TU C VU* C ), H {-aj). 
This obviously explains why the overall sign of the spectral shift function is reversed under the change A 1— > 
—A. But it also explains why the roles ofVu an d V33 are interchanged in the estimates. Indeed, the natural 
projection corresponding to the vector potential a is P = diag(P, 0, P. 0) since we have bo > Ofor a, whereas 
P' := diag(0, P, 0, P) is the natural choice for the vector potential —a since we have bo < Ofor —a. Now, one 
has 

V44 — V43 — V42 V41 
~ V34 V33 V32 — V31 
— V24 V23 V22 —V21 

V14 -Vii"-Vi 2 V11 

So, f/je projection P which selects the components ±(Vii,V33) of the potential ±V is replaced, after the 
change A 1— > —A, fry f^c projection P' which selects the components ^(^33) Vn) of the transformed poten- 
tial zpUcVUc. 

For the next proposition we define for each A G Kwith|A| > m the positive operator fiW (A) in L 2 (R 2 ; C 8 ) 
given by 



fi^(A) := —7=!== fP^ xp °) wh ere M, 



(|A+to|W+ 0\ 
1 
o' A "™'^o/' 

Proposition 6.16. (a) Let V satisfy Assumption 3.2 with v £ (3, 4]. Then one has for each s > and each 

TV { arctan [s _1 f2(A)] - arctan [s" 1 ^ 1 ' (A)] } = 0(| A T m\~ 5 ) as A — > ±m, |A| > m. 

(») Lef V satisfy Assumption 3.1 with v± > 2 anrf 1/3 > 2. 77ien one has for each s > 

Tr{ arctan [s" 1 fi(A)] - arctan [s _1 Q (1) (A)] } = 0(1) as A -> ±m, |A| > m. (6.27) 

Proof. Points (a) and (b) are proved by using the Lifshits-Krein trace formula (1.2) with /(A) = arctan(A), 
AeM. We do not give the details, since the argument is analogous to the one of [12, Cor. 2.2]. □ 

Note that if V satisfy Assumption 3.2 with v £ (3, 4], we can choose 5 <G (^^j ^zr)> an d so Proposition 
6. 16. (a) entails 

Tr{ arctan [s -1 fi(A)] - arctan [s" 1 ^ 1 ' (A)] } = o(|A =p m|"^i) as A — * ±m, |A| > m. (6.28) 

Moreover, if V satisfy Assumption 3.2 with 1/ > 4, then it satisfies Assumption 3.1 with u± > 2 and 1/3 > 2, 
and, hence (6.27) is valid. Finally, we have for s > and |A| > m 

Tr arctan (A)] = jf^ n + (2st(^) 1/2 ; P W + p) + jf °° ^ n + ( 2s ^) 1/2 ;^ 

(6.29) 

Combining Equations (6.27)-(6.29), Theorem 6.14 and Lemmas 6.7-6.9, we get the following. 
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Corollary 6.17. Let V satisfy Assumption 3.2. 

(a) Assume that the hypotheses of Lemma 6.7 hold with U± — W± and a = v — 1. Then we have 

C(A; g_, g ) = - - cqs ^ _ 1)} (2(^) 1/2 ; u + , b ) (l + o(l)) as A \ m, 

£(A; #+,#„) = - ^ — ^_ 1 (2(^) 1/2 ; U _,6o)(l + (l)) as A / -m, 

2 cos (tt/ (i/ — ljj V / 

vv/f/i Vt^-i given by Equation (6.21). 

(b) Suppose that V also satisfies (3.1) with v±_ > 2 a«t/ 1/3 > 2, and assume that the hypotheses of Lemma 

6.8 hold with U± = W±. Then we have 

C(A;g_,g ) = -^/ 3+ (2(^) 1/2 ;77+^o)(l + o(l)) as A \ m, 

^(A;g_,g ) = i$ /3 _(2(^±A) 1/2 ; ?? _,& ) (l + o(l)) as A / -m, 
vw'f/i /3± € (0,oo) and&/3 ± given by Equation (6.22). 

(cj Suppose that V also satisfies (3.1) vv/f/z > 2 a«c/ 1/3 > 2, a«t/ assume that the hypotheses of Lemma 

6.9 hold with U± = W±. Then we have 

e(A;g_,g ) = -i<i> 00 (2(^) 1/2 )(l + (l)) as A \ m, 

^(A;g + ,g ) = i$ oo ( 2 (i|±A) 1/2 )(l + o (l)) as A / -m, 

vv/f/z $oo given by Equation (6.23). 

Putting together the results of Proposition 6.10 and Corollary 6. 17, we obtain the following. 

Corollary 6.18. Under the assumptions of Corollary 6. 17. (a), we have 

nm Z(mQ.-e)- 1 ;H-,Ho) = 1 = g( - m(l - e)' 1 ; H+, Hp) 

e\o £(m(l-e);g_,g ) 2cos (tt/O - 1)) e\o £(-m(l- e);H+,H ) ' 

and under the assumptions of Corollary 6.17.(b)-(c), we have 

CHl-e)- 1 ;^,^) = 1 = ^-mll-e)-';^,^) 
e™ eMl-£);g-,g ) 2 e™ e(-m(l-£);g+,g ) 
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7 Appendix 

We give in this appendix the proof of the inclusion (3.6) for the class of potentials V given in Remark 3.3. We 
start with a technical lemma. We use the notations a := (ai, 02, a 3 ) T and 

(d t v) := {(d t v jk )}, w := (fliV.ftv; W", (did m v) ~ {{d,d m v ok )}. 

Lemma 7.1. Let V be as in Remark 3.3. Then 
(a) One has in 3§(T>(Ho), 2?(i?o)*) the equalities 

[H ,H] = -ia • (W) + [03, y]P 3 + m[/3, V] (7.1) 
= -ia • (VI/) + P 3 [a 3 , V] + i[«3, (9 3 y)] + m[/3, V]. (7.2) 

(£>) Lef z € R \ {cr(-ffo) U a(H±)}. Then there exist operators B± G £§(7i) such that 

R 2 ±{z) = B ± Hq 2 and R\(z) = Hq 2 B* ± . (7.3) 

Proof, (a) We know from Lemma 2.2(b) of [32] that V(Hq) C X>(P 3 ). So each member of Equations (7. l)-(7. 2) 
belongs to 3§(V(H ), V{H Q f). 

Let ip G V(Hq), take a sequence {</2„} C Cg^M 3 ; C 4 ) such that lira ||y? n — (fi\\v<H ) = 0, and denote by 

n 

(-, the anti-duality map between V(Hq) and V(Hq)*. Then 

(vj, [fr 0j f?]v»)i,-i = ^) - {v<p,H (p) 

= lim. (ip n , [ai(Pi - ai) + a 2 (P2 - 02) + a 3 P 3 + /3m, V]<p n ) 

= lim(p„, { -ia- (W) + [a 3 , V]P 3 + m[/3, V]}cp n ). (7.4) 

Since T>{Hq) C 2?(P 3 ), we also have lim \\ip n — <p\\t>(p-,) = 0, and thus 

n 

(p, [#o,#Mi,-i = (p, { - ia ■ (VF) + [a 3j ^]P 3 + m[/3, V]}ip) 

This proves (7.1). Using (7.4), one also gets the equality (7.2). 

(b) In what follows, we omit the indices "±" to simplify the notations and we write B±, B2, ■ ■ . for elements 
of 3s\hL). Since V(H) = V(H ), we have 

R 2 {z) = BiHq 1 R(z) = BiRiz)^ 1 +B 1 [Ho\R(z)] = B 2 H« 2 + B^ 1 R(z)[H , H]R{z)H^ . 

Now, one has 

R{z)[H Q ,H]R(z)H^ = R(z){-i a -(VV)+P 3 [a 3 ,V]+i[a 3 ,(d 3 V)]+m[(3,V]}R(z)H^ = B 3 H„ 2 

due to Equation (7.2), the equality T>(H) = T)(Hq), and the inclusion T>(Hq) C T>{P 3 ). This, together with 
the preceding equation, implies the first identity in (7.3). The second identity follows from the first one by 
adjunction. □ 

Proposition 7.2. Take z€t \ {o-(Hq) U <j(H±)} and let V be as in Remark 3.3. Then we have 

R 3 ± {z) - R 3 (z) G Si(H). 

Proof. In what follows, we omit the indices "±" to simplify the notations and we write B\ , P2 , ■ ■ ■ for elements 
of £$(TL). Differentiating twice the resolvent identity 

R(z) - Ro(z) = -R{z)VRq{z) 
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we find that 

R 3 {z) - Rl(z) = -R{z)VRl{z) - R 2 {z)VRl(z) - R 3 {z)VR {z). 

So it is sufficient to show that each term on the r.h.s. belongs to Si(TC). This is done in points (i), (ii) and (iii) 
below. 

(i) For the term R(z)VRf ) (z), one has 

R(z)VR 3 (z) = R(z)R (z)VR 2 (z) + R(z)[V, R (z)]R 2 (z). (7.5) 

Since V(H) = V(H Q ), one has 

R(z)R Q (z)VR 2 (z) = R(z)(H - z)R 2 (z)V R 2 {z) = (B^ 2 V 1 ' 2 ) (V 1 ' 2 H^ 2 B 2 ) . 

So, by (3.3), R(z)Ro(z)VRq(z) is the product of two Hilbert-Schmidt operators, and thus belongs to Si(H). 
For the second term of (7.5), we have by (7.1) 

R(z)[V,Ro(z)]R 2 (z) = R(z)R Q (z)[H Q ,H]R 3 (z) = B 1 H- 2 { - ia • (W) + [a 3 ,V]P 3 +m[P,V]}H^ 3 B 3 . 

Due to the hypotheses on V and (djV), one can use (3.3) to write the first and third term as a product of two 
Hilbert-Schmidt operators. So it only remains to show that H^ 2 [a 3 , V]P 3 Hq 3 belongs to Si(H). For this, we 
use the inclusion T>(Hq) C T>(P 3 ) and the commutation of P 3 and Hq 1 on T>(P 3 ) [32, Lemma 2.2(b)] to get 

H^ 2 [a 3 ,V]P 3 H^ = H^ 2 [a 3 ,V]H^ 2 P 3 H^ = H^ 2 [a 3 ,V]H^ 2 B A . 

This, together with (3.3), implies that H^ 2 [a 3 ,V]P 3 H^ 3 belongs to Si(H). 

(ii) One can write R 2 {z)VR 2 (z) as the product of two Hilbert-Schmidt operators by using (7.3) and (3.3): 

R 2 {z)VR 2 {z) = B 5 H^ 2 VH^ 2 B 2 = (B 5 H^ 2 V 1 ' 2 )(V 1 ' 2 H^ 2 B 2 ). 

Thus R 2 {z)VRl{z) belongs to S^H). 

(iii) For the term R 3 (z)V Rq(z) we have 

R 3 (z)VR a (z) = R 2 (z)VR(z)R (z) + R 2 (z)[R(z),V}R (z). 

One shows that R 2 (z)V R(z)R^(z) € Si(H) as in point (ii). For the second term, we have by (7.2) and (7.3) 

R 2 (z)[R(z),V]R (z) = R 3 (z)[H,H Q ]R(z)R Q (z) 

= B 5 Hv 2 R(z){ia ■ (W) - P 3 [a 3l V] - i[a 3 , (d 3 V)} - m[(3, V]}H„ 2 B 6 . 

Due to the hypotheses on V and (djV), one can use (7.3) and (3.3) to write the first, third, and fourth term as a 
product of two Hilbert-Schmidt operators. So it only remains to show that H^ 2 R(z)P 3 [a 3 , V]Hq 2 belongs to 
Si(H). Using [32, Lemma 2.2(b)] and (7.3), one gets 

H^ 2 R(z)P 3 [a 3 , V]H^ 2 = H^ 2 P 3 R{z)[a 3l V]H^ 2 + H^ 2 [R(z), P 3 ][a 3l V]H^ 2 

= P 3 H^ 2 R(z)[a 3 ,V]H^ 2 - iH^ 2 R(z)(d 3 V)R(z)[a 3 ,V]H^ 2 
= B 7 H^ 2 [a 3 ,V]H^ 2 + B s R(z)(d 3 V)R(z)[a 3 ,V]H^ 2 . 

The first term on the r.h.s. belongs to Si(H), and for the second term we have by (7.2) and (7.3) 

R{z){d 3 V)R{z)[a 3 ,V]H« 2 

= (d 3 V)R 2 (z)[a 3 ,V]H^ 2 + R(z)[(d 3 V),H}R 2 (z)[a 3 ,V}H^ 2 
= B 9 H^ 2 [a 3 , V]H^ 2 + R{z){ia ■ [V(d 3 V)\ - P 3 [a 3 , (d 3 V)} 

-i[a 3 ,(d 2 V)]-m[p,(d 3 V)} + [{d 3 V),V]}B b H^ 2 [a 3 ,V]Hv 2 . 

Due to the hypotheses on V, (djV), and (dj 3 V), one can use (3.3) to show that the first, the second, the fourth, 
the fifth, and the sixth term are trace class. For the third term we have to use (3.3) and the fact that R(z)P 3 
extends to a bounded operator. □ 
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Remark 7.3. When the potential V is scalar, the equations (7.1)-(7.2) reduce to the single equality 

[H ,H) = -ia-{VV) 

in 3$(T>(Hq), D(H )*). So the calculations in points (i) and (Hi) of the proof of Proposition 7.2 simplify ac- 
cordingly, and we obtain the inclusion 

R\{z) - Rl{z) e Si{H) 
without assuming anything on the derivatives ofV of order 2. 
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